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We propose a new method to evaluate the Lellouch-Liischer factor which relates 
the AI = 3/2 K — > tttt matrix elements computed on a finite lattice to the physical 
(infinite- volume) decay amplitudes. The method relies on the use of partially twisted 
• boundary conditions, which allow the s-wave tttt phase shift to be computed as an 
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almost continuous function of the centre-of-mass relative momentum and hence for 
its derivative to be evaluated. We successfully demonstrate the feasibility of the 
technique in an exploratory computation. 

PACS numbers: 11.15.Ha, 12.38. Gc, 13.25.Es 

I. INTRODUCTION 



X 

The precise computation of K — > txtx decay amplitudes in lattice simulations would be 



a very important step in our ability to evaluate non-perturbative strong interaction effects 
in hadronic processes. In particular it would allow us to compute the quantity e'/e which 
contains the QCD effects in direct CP-violating decays and to understand the origin of the 
long-standing puzzle of the AI = 1/2 rule. Such calculations have not yet been achieved, 
primarily due to the technical difficulties in computing the AI = 1/2 K — > tttt matrix 
elements and the corresponding disconnected diagrams with sufficient precision. 

In the centre-of-mass frame, the relation between K — > tttt matrix elements computed 
in a finite Euclidean volume and the physical decay amplitudes is given by the Lellouch- 
Liischer (LL) factor . This relation was generalised to a moving frame, in which the 



total momentum of the two-pion system is non-zero, in The LL factor depends on 

S'(q*), the derivative of the s-wave phase-shift S(q*) with respect to the relative momentum 
in the centre-of-mass frame q*. Whereas the phase-shift itself can be determined from the 

nrr 

two-pion spectrum in the finite volume p, |6[, the derivative can only be determined by an 
interpolation or from estimates using chiral perturbation theory or models. The aim of this 
paper is to demonstrate that for AI = 3/2 transitions it is possible to calculate both the 
phase-shift and its derivative directly at the quark masses used in the lattice simulation. 
This can be achieved by using (partially) twisted boundary conditions. 

This paper is organized as follows. In the next section, we summarise the theoretical 
background, recalling the expression for the LL-factor and describing the main properties 
of partially twisted boundary conditions. In this section we also present our procedure for 
the determination of the derivative of the two-pion phase shift. We then test our proposed 
procedure in a numerical simulation which is described in sec JIIII and the results are presented 
in sec JIVI Finally in sec]V]we present our conclusions. 

II. THEORETICAL FRAMEWORK 

We start with a reminder of why it is necessary to compute the derivative of the phase 
shift. In an arbitrary moving frame the K t\t\ matrix elements in finite and infinite 
volumes, which we denote by M and A respectively, are related by jl-4] 

\A\* = 8nV*^f (W)+0 P V)} \M\\ (1) 

where E is the total energy, E 2 — P 2 = 4(m 2 + q* 2 ) where P is the total momentum, the 
' denotes differentiation with respect to q*, 5 is the s-wave irir phase-shift and V = L 3 is 
the spatial volume. The explicit form of the kinematic function <p p can be found in 4] 
and we do not need to reproduce it here. Our aim in this paper is to explain how S'(q*) 
can be calculated directly for AI = 3/2 transitions and to demonstrate the feasibility of the 
proposed procedure in an exploratory lattice computation. In this way the factor relating 
the finite and infinite volume matrix elements can be determined at the quark masses used 
in the simulation. 

In order to evaluate the 1 = 2 s-wave phase-shift and its derivative we compute the cor- 
relation function {0\7T~ (t)n~ (t) 7r + (0)7r + (0)|0), where ir + and ir~ are interpolating operators 
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with the correct quantum numbers to create or destroy a tt + meson. The time t is chosen to 
be positive. The precise form of the interpolating operators used in our exploratory simula- 
tion will be described in sec. IHII below; the general argument presented here does not depend 
on the choice of these operators. With periodic boundary conditions the total momentum of 
the two-pion system P can take the values P = (2ir/L) n, where n is an vector of integers. 
From the energy levels of the two-pion system we can then deduce the corresponding phase 
7| using Liischer's quantisation formula: 



shifts 



*(?*) + <lf{<f) = rwr , (2) 

where n is an integer. In typical simulations 2n/L is several hundred MeV and so the allowed 
momenta are quantized in large increments making it impossible to calculate reliably the 
derivative of the phase-shift. 

Changing the boundary conditions changes the momentum spectrum and Bedaque pro- 
posed exploiting this fact to extend the range of momenta available in lattice applica- 
tions 8j 1 . For example if we impose the following spatial boundary condition for the 
w-quark: 

u(xi + L) = e i6i u{xi) , (i = l,2,3) (3) 
then the momentum spectrum of a free ir + meson is given by 

27r , 9 i fA^ 

Vi = rii — + — , (4) 

where the rij are integers. For illustration let us imagine that only 9\ ^ and we now drop the 
suffix and denote 9\ by 9 (the generalization to arbitrary 9 is conceptually straightforward). 
For 9 < ir the ground state of the two-7r + system then has momentum 29/ L, corresponding 
to each pion having momentum 9/L. Boosting to the centre-of-mass frame we find that 
q* = for all 9 so that nothing has been gained by introducing the twisted boundary 
conditions. This is not true if we add some units of 27r/L to the momentum by performing 
the corresponding Fourier transform. In this case, q* does depend on 9 so that by measuring 
the two-pion energies with different 9 we obtain the phase-shifts at different values of q*. 
Since 9 can be incremented by arbitrarily small amounts, in this way we can obtain the 
derivative of the phase shift. This is the main point which we wish to make in this paper. 



1 See [8( and [9( for references to earlier papers with related ideas. 
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As an illustrative example let the w-quarks satisfy twisted boundary conditions with 
angle 9, the d quarks satisfy periodic boundary conditions and the two-pion state be given 
a further momentum of — 2tt/L by performing the corresponding Fourier transform. The 
total momentum of the two-pion state is then Pg = 2(9 — ir)/L. We imagine that the matrix 
elements of the operators appearing in the AI = 3/2 weak Hamiltonian have been measured 
with periodic boundary conditions (i.e. with 9 = 0) with a total momentum of —2n/L and 
that we want to evaluate the corresponding LL factor. We now measure the two-pion ground 
state energies Eg for a range of values of 9 and determine the corresponding centre-of-mass 
relative momentum q* in the standard way using Eg — Pg = 4(m^ + Qg 2 )- The phase shift 
S(qg) is obtained from the generalization of the Liischer quantization condition to a moving 
frame ((2]) js], 0, [lO] and its derivative is determined from the observed slope of the results 
obtained at different values of 9 which can be chosen to be arbitrarily close together. 



The met 
conditions 



rod proposed in this paper is feasible with the use of partially twisted boundary 



11 



12j in which the twisted boundary conditions are applied only to the va- 



ence quarks, whereas the sea quarks satisfy periodic boundary conditions. It was shown in 



11] that with partially twisted boundary conditions the finite-volume effects remain expo- 



nentially small. The practical advantage of not varying the boundary conditions for the sea 
quarks is that it is not necessary to generate a new ensemble of gauge configurations for every 
choice of twisting angle; indeed it is this feature which makes the method feasible. Partially 
twisted boundary conditions have been shown to satisfy the dispersion relation and used in 
the evaluation of the leptonic decay constant f n 13| , in the evaluation of the electromagnetic 
form factor of the pion and of K — > n semileptonic decay amplitudes [171 ]. 

It is convenient to implement the partially twisted boundary conditions, by introducing 
a change of quark field variables ip — > 



ip(x) = e l ~ip(x) and ip(x) 



'tp(x) 



(5) 



so that the ip(x) satisfy periodic boundary conditions. The hopping terms in the lattice 
fermion action now become (for i = 1, 2, 3) 



ii>(x) e tZ ?Ui{x){l -7i)V>(x + z) + e ^U}(x - i)(l + 7^0 - «) 



(6) 



The quark propagator is obtained by inverting the Dirac operator in the gauge field back- 
ground {e l ^Ui(x)} where the {Ui(x)} are generated from a dynamical Markov chain with 
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sea quarks. The resulting Dirac operator is hermitian and standard conjugate gradient 
algorithms can be used to invert it. 

We end this section with three simple observations. 

i) With the prescription described above the phase shifts obtained at different values 
of q* are correlated which, as explained in section IIVBI allows for a more precise 
determination of the slope than would be the case if the errors on the individual 
points were independent of each other. 

ii) The second observation is that the choice of the state is the most convenient 
one with which to obtain 1 = 2 phase-shifts, in spite of the fact that kaons do not 
decay into 7r + 7r + . The tx + ti + state automatically has 1 = 2 whereas, because the 
twisting cancels in the 7r°, it is not possible to symmetrize the ir + 7T° state to separate 
the 1 = 2 and 1=1 components. The Wigner-Eckart theorem allows us also to obtain 
K + — » 7r + 7r° matrix elements of the operators in the effective weak Hamiltonian by 
computing K + — > 7i + tt + matrix elements of the corresponding AI Z = 3/2 operator: 

4 {<T + (PiK(fc) I + (n + {p 2 )APi) 1} |0 3/2 | K + ) = l{n + ( Pl )n + (p 2 ) |0' 3 / 2 | K + ) . 

(7) 



O 3 / 2 is one of the three A J = 3/2 operators: 



= (?<T)l {(W)l - {d*d?) L } + (5V)i {v?d?) L (8) 
OT = (s l d l ) L {(u j u j ) R - (d j d j ) R } + (sV)x (u j di) R (9) 
Ol' 2 = (s\J0l {(«''«%- [d^) R } + {tu^L^d^R, (10) 

where i and j are colour indices, {q~\q2)L,R = 9i7 M (l T 7 5 )?2 and O' 3 ^ 2 is the corre- 
sponding operator with AI Z = 3/2: 

'(2 3 7 /2 i) = )^ ^ dj )L , 0' 7 3/2 = (s l d*) L (uidJ) R , 0' 8 3/2 = (s l di) L (u j d% . (11) 

iii) Finally we recall that a finite cubic box breaks rotational invariance and hence mixes 
different partial waves. For lattice 7r + 7r + states with partially twisted boundary con- 
ditions in addition to the s-wave, there are also components with 1 = 2,4,... and the 
results presented here were obtained under the assumption that only scattering in the 
s-wave is important in the energy range of interest (which is a good approximation for 
K — > hit decays). 
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This concludes the demonstration that, in principle at least, it is possible to calculate 
the derivative of the s-wave nn phase-shift, and hence the Lellouch-Liischer factor, at the 
masses and momenta used to calculate the finite- volume K — > titt matrix element. In the 
following section we test this idea in an exploratory simulation. 



III. DETAILS OF THE LATTICE COMPUTATION 



We now report on our exploratory numerical study of the ideas proposed in the previous 
section. The computations are performed on a 16? x 32 x 8 lattice with Nf = 2 + 1 Domain 



Wall Fermions and the Iwasaki gauge action [18] a t = 2.13. Details of the simulation 
and properties of the ensembles are described in 19J. The bare strange quark mass is fixed 
at am s = 0.04 and we use two ensembles with light quark masses am u = am^ = 0.04 and 
0.02. (We will denote the common light-quark mass by m and the strange quark mass by 
m s .) The lattice spacing was estimated in 19[ from the static potential using the value 
ro = 0.5fm and was found to be 1/a = 1.826(48)(94) GeV (in the chiral limit). Since our 
goal is to demonstrate the feasibility of the method rather than to obtain precise physical 
results we will use this value whenever we want to convert our results to physical units. We 
find that the "pion" masses are am n = 0.364(2) for am = 0.02 and 0.468(2) for am = 0.04 2 . 



There are 2595 thermalized trajectories generated for each ensemble [19J]. For this study 
when m = 0.04 we use 250 configurations separated by 10 trajectories, whereas for m = 0.02 
we use 120 configurations separated by 20 trajectories. In order to reduce the statistical 
error, for each configuration we measure correlation functions with two different sources (at 
t = and t = 16). In addition for each choice of 9 = \9 | we perform measurements on 
each configuration with 9 = (0,0,0), (0, 9, 0) and (0,0,0). The autocorrelations are studied 
using the standard binning technique, in which the results from a number of neighbouring 
configurations are combined in bins and treated as independent measurements. The bin size 
is increased until the corresponding statistical error stabilizes; for our observables this is the 
case with bins of size 4 for m = 0.02 and 7 for m = 0.004 (recall that the configurations 
are separated by 20 and 10 trajectories respectively for the two masses). The bin sizes are 
a little larger than those found for the observables studied in 19] . 



2 In order to simplify the notation, in the remainder of the discussion we use lattice units setting a = 1. 
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A. Interpolating operators 

For the interpolating operator for each pion at the source we use the gauge-fixed wall 
source with momentum projection, 

Ol(p u ,p d ) = ^u@)A\2)^* 7s '£ i Atf)dtf)t&*, ( 12 ) 

x y 

where A(x) is the gauge fixing matrix. In this study we use the Coulomb gauge. The 
superscript i in 0\ stands for initial. At the sink we use a point source 

Ol(p) = ^u{x) l5 d{xy?- s . (13) 

x 

The two-pion interpolating operators are then constructed as the product of the single pion 
ones, 

0^ = 0^ u ,p d )0^ u ,p d + P) and Ol = Oi(p)Oi(p + P) (14) 

— * 

where P = (2n/L) n and n is a vector of integers. We choose both pions to be created (or 
annihilated) on the same time slice. 

B. Twisted boundary condition 

In this calculation partially twisted boundary conditions with twisting angle 9 are applied 
to the w-quarks, while the rf-quarks satisfy periodic boundary conditions. Thus, the momenta 
of free 7r + -mesons are quantized as 

9 2tt^ 

p w + = — + — n . (15) 

We choose to introduce momenta only in a single direction, the ^-direction say, and take a 
fixed P = (0,0, -2-k/L) for a variety of values of 9 (9 = (0,0,9)) with < 9 < ir/L. The 
choice of | P\ to be the smallest non-trivial Fourier momentum is motivated by the need to 
keep lattice artefacts as small as possible. 

An advantage of the change of variables introduced in Eq. §5§ is that the phase factor 
for momentum 9/L is already absorbed in the quark field. Therefore, the same form of 
interpolating operator can be used for all twisting angles. 
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(c) (d) 

FIG. 1: The four diagrams contributing to the ir + ir + correlation function. The grey ovals represent 
the gauge-fixed wall sources for each of the ir + interpolating operators and the black circles represent 
the local interpolating operators at the sink. 

C. The Diagrams 

In the construction of the correlation function (0\OH(t) 0^.(0) |0), where only the time 
variable is explicitly shown, there are four diagrams (or Wick contractions) as indicated in 
FigfTJ The two disconnected diagrams (D) in the first row of FigJT] each have two traces 
over the spinor and colour labels and for each gauge configuration are the product of two 
separate pion loops. The two connected diagrams (C) in the second row of the figure have 
a single trace over the spinor and colour labels and so the two-pion correlation function is 
the difference D — C, where the — sign comes from Fermion statistics. 

IV. RESULTS 

In this section we present the results of our numerical study. We measure the energies 
of the two-pion states for a range of values of 9 and then use the Liischer quantization 
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FIG. 2: Dispersion relation for pseudoscalar mesons with light quark masses m u = 0.02 and 0.04. 



condition to determine the corresponding 7T7r-phase shift. By choosing the twisting angles 
to be sufficiently close together we are then able to calculate the derivative of the phase 
shift and hence the corresponding LL factor. We start however, with a demonstration of the 
dispersion relation for a single pion as a function of the momentum induced by the twisting 
angle, p = 8/L. In Figj2] we plot the measured value of E 2 as a function of p 2 = \p\ 2 
and compare it to the continuum E 2 = p 2 + m 2 dispersion relation as well as a possible 



discrete form, E 2 



m 



sin 2 (jo). The mass of the pion is determined in the standard way 



by fitting the correlation function at zero momentum so that the theoretical curves have no 
free parameters. The results are presented for two different values of the bare light quark 
mass m = 0.02 and 0.04 corresponding to pion masses m n = 0.364(1) and m n = 0.4676(8) 
respectively (to convert the pion msses into physical units, recall that the lattice spacing a 
obtained from the Sommer scale is approximately given by 1/a = 1.83(5)(9)). In ref. 19] 
the residual mass for this action was found to be approximately m res ~ 0.011 and using this 
value the ratio of pion masses is as expected from the PCAC relation. 

As is now well established (see for example ref. 13[), the dispersion relation is well sat- 
isfied. The small difference in the theoretical curves at larger values of the momentum is 
due to 0((ap) 4 ) terms. As the momenta approach p 2 = (2n/L) 2 ~ 0.154, the smallest mo- 
mentum which is accessible with periodic boundary conditions, there is a small discrepancy 
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visible between the results of the lattice calculation and the continuum dispersion relation 
for the m = 0.04 case. For the purposes of this feasibility study we neglect this difference. 

A. The energies of the two-pion states and the corresponding phase shifts. 

In this exploratory study we calculate the derivative of the phase-shift by using the simple 
ansatz 

dS_ ^ AS_ = 6(4) - SjgU) 

dq* Aq* q* - q\_ x ' ^ ' 

where q* and q*_ l are neighbouring values of q*. Since both the numerator and denominator 
are small quantities, it is necessary to exploit the fact that all the quantities on the right- 
hand side of ffTB"]) are calculated on the same configurations and are therefore correlated. 
This can be achieved by fitting all the necessary quantities, i.e. the q* and corresponding 
phase-shift, simultaneously. 

We start by considering the interaction energy 

AE = E T7T - y/ml + |pi| 2 - yjm\ + |p 2 | 2 , (17) 

where E nn is the energy of the two-pion state and p\ and p 2 are the corresponding momenta 
of each of the pions if there were no interactions between them. The results for AE for 
m = 0.02 obtained by measuring E nn and m n and boosting to the centre-of-mass frame are 
shown as the red squares in FigfSJ From our numerical study however, we find that the 
errors in AE can be significantly reduced if instead of using eq. (IT7|) we use 

AE = E n7T — E n> i — E n> 2 , (18) 

where E nA and E v ^ 2 are the measured energies of free pions with momenta pi and p 2 respec- 
tively. By performing a jackknife analysis of the right- hand- side of ( TlBl) . we can take the 
correlations between E nn and E^i + E n ^ into account. The results are plotted as the black 
circles in FigSJ, from which we see that the errors are decreased by more than a factor of 2. 
As the input two-pion energies into the Liischer quantisation formula we therefore take 

EZ = {E^ - K,i - K l2 } + V^4 + W+ + IP2I 2 . (19) 
The single-pion correlation functions are fit to the standard form 

C^(t) = B cosh(£^(t - T/2)) , (20) 
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FIG. 3: The energy shift in the centre-of-mass frame as a function of q* obtained using eq. (|17|) (red 
squares) and (|18p (black circles). The graph shows the advantage of exploiting the correlations in 
the two-pion and single-pion correlation functions. 

where % = 1 or 2 and T = 32 is the temporal extent of the lattice. 

The leading behaviour of the two-pion correlation function takes the form: 

C^{t) = Acosh(E nn (t - T/2)) + Ccosh(D(t - T/2)), (21) 

where the second term on the right-hand side is the contribution in which each of the two-pion 
interpolating operators destroys one pion and creates another, so that one pion propagates 
from to t and the other from t to T. Denoting the energies of the two pions in the second 
term by E Wi x and then D = \E njl — E ni2 \ and C is proportional to e -{ E ^^+ E ^a) T / 2 . j n 
the first term A is proportional to e~ E7TwT ^ 2 and so, except in the vicinity of t = T/2, the 
first term in (l2Tj) dominates the second. 

By simultaneously fitting C^i(t), C^^if) and C W7T (t) to the forms in eqs. ( 120]) and ( 12T|) all 
the necessary quantities can be determined. We find however that the errors are reduced 
if we simultaneously fit the correlation functions Cn,i(t), C^pit) and the ratio i? W7r (t) = 
C nn (t) / C ni i(t)C ni 2(t) and the results quoted in the tables were obtained in this way. (To 
speed up the minimization of x 2 ; we found it useful to first fit C n> i(t), C^^if) and CW^) 
and then to use the resulting parameters as the starting values in the fits for C nt i(i), C n ^(t) 
and Rn n (t).) 



12 



a 

V 


PL 




A TP n* 

L\Hj q 


ri=2 


n 
u 


2vr 0.9220(36) 


u.uziu^iuj u.zuiy^izj 


n 97/1 fi ^ 


7T 

18 


17 TV 

o 
y 


0.9145(36) 


U.UZlo^lUJ U.ZUOf^llJ 


n 977f 1 9 s ! 
-U.Z 1 f ^1Z ) 


7T 

9 


16tt 
o 


0.9075(35) 


n n9i fii fQ^\ n9n c ;^i'in s i 


n 98iV 1 1 ^ 

-U.ZoU^ll ) 


7T 

6 


57T 

Q 
O 


0.9008(35) 


n 091 QltQl \ D9n79i'in s i 
u.uziyo^yij u.zuiz^iuj 


n 989f1 1 
-U.ZoZ^ll J 


2- 

9 


14?r 

o 
y 


0.8946(34) 


nn999/i/'87\ n 9nssoi'Q c i'i 
u.uzzz^o ( i u.zuooy^yoj 


n 98^^ 1 n\ 


7T 

4 


3tt 

n 
z. 


0.8917(34) 


u.uzzoy^ooj u.zuyoy^yzj 


n 987^ 1 

-U.Zo f I^IU ) 


5- 
18 


13tt 

Q 

y 


0.8888(34) 


099^/1 n 91 f\A7(RQ\ 
U-UZZO^ooJ U.ZlU4t(^oyj 


n 988^ 1 

-U.Zoo^lU J 


7T 

3 


i7T 

3 


8835(34) 


n n99S' : i('8n N \ n 91 1 Qfii's c i s i 
u.uzzoo^ouj u.ziiyo^ooj 


n 9QD7l'Q7 

-u.zyu / ^y / 


7T 

2 


7T 


0.8700(35) 


U.UZoOo^/OJ U.ZlOoO^/OJ 


n 9Q7Ql'8Q 

-u.zy * y ^oy 


2- 
~3~ 


2tt 
3 


0.8604(33) 


0.02415(71) 0.21874(69) 


-0.3035(84 


3tt 
4 


7T 

2 


0.8570(32) 


0.02435(71) 0.21977(68) 


-0.3055(83 


5- 
6 


7T 

3 


0.8546(29) 


0.02450(73) 0.22052(68) 


-0.3070(84 


7T 





0.8526(29) 


0.02465(76) 0.22115(71) 


-0.3086(88 



TABLE I: Results for the total two-pion energy energy shift AE, the relative centre-of- 

mass momentum q* and the 1 = 2 s-wave phase shift <5j =2 for quark masses m/m s = 0.02/0.04. 
The total momentum P = 2(ir — 9)/L. 

The fit ranges for C % i(t) and C^^if) have already been determined from the standard 
effective mass plots. For the two-pion correlation functions C n7T (t) we find that the range 
t G (8, 17) is consistent for all twisting angles and for both values of the quark mass and we 
use the same range for R^(t). Since our final goal is to determine 85 /dq*, having a single 
fitting range for all the twisting angles simplifies the analysis. 

The energies of the irir states are given in Tablesjl] and [III for the two values of the light 
quark mass. In these tables we also give the values of the energy shift AE, the relative 
momentum in the centre-of-mass frame q* and the phase shift. The state with the largest 
total momentum has the largest total energy as expected, but it also has the smallest relative 
momentum q* because of the boost factor. In FigJH and Fig[5j we plot AE in the centre-of- 
mass frame and the phase-shift as a function of q*. AEqm is obtained from the measured 
value of AE in the moving frame (see eq. (118j) ) by AEqm = 1 AE, where the boost factor 
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1.0351(29) 


0.01777(42) 


0.21754(51) 


-0.2736(62) 


5tt 
6 


7T 

3 


1.0332(29) 


0.01788(41) 


0.21808(49) 


-0.2749(61) 


7T 





1.0318(29) 


0.01803(40) 


0.21859(48) 


-0.2769(59) 



TABLE II: Results for the total two-pion energy energy shift AE, the relative centre-of- 

mass momentum q* and the 1 = 2 s-wave phase shift <5j =2 for quark masses m/m s = 0.04/0.04. 
The total momentum P = 2(ir — 9)/L. 



7 = 1 / Vl — v 2 and v = \pi +p 2 \/E 7Tn . As expected, the energy shift E C m is positive and 
increases as the relative momentum q* increases. The phase shift is negative reflecting the 
repulsive nature of the 1 = 2 two-pion interaction. 

B. The derivative of the phase shift and the Lellouch-Liischer Factor 

Once the phase shift has been determined as a function of q*, the derivative can be 
calculated and the Lellouch-Liischer factor can be estimated. The values of the phase shift 
at different q* are correlated of course and, as we shall see, this results in a smaller error 
on the slope than might have been expected simply by looking at the error bars on a 
pair of neighbouring points and considering them to be independent. This is because the 
results for the phase-shifts at different q* move up and down together as we sample different 
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FIG. 4: The energy shift in the centre-of-mass frame as a function of q* at the two values of the 
light-quark masses. 
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configurations, leading to a smaller variation of the slope. 

There are many ways in which the derivative can be determined numerically. We choose 
to use the most simple-minded one, i.e. the simple ansatz in eq. (fT6]) . Note that the ansatz 
provides the derivative most accurately at the mid-point between and qi-\. We therefore 
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q* d5/dq* d(j) p /dq* LL factor 



0.2028(11) 


-1.49(54) 10.873(13) 4027(224) 


0.2046(11) 


-1.55(46) 11.000(12) 3928(189) 


0.2064(10) 


-1.62(45) 11.126(11) 3835(179) 


0.20807(98) 


-1.68(43) 11.250(10) 3751(172) 


0.20929(93) 


-1.73(42) 11.3707(92) 3704(166) 


0.21008(91) 


-1.75(42) 11.4297(87) 3669(165) 


0.21122(87) 


-1.79(42) 11.4867(82) 3609(162) 


0.21391(80) 


-1.85(40) 11.5953(73) 3472(153) 


0.21730(72) 


-1.93(47) 11.8962(54) 3368(176) 


0.21926(68) 


-1.96(63) 12.1328(46) 3335(227) 


0.22015(68) 


-2.02(83) 12.2188(45) 3296(294) 


0.22084(70) 


-2.4(1.5) 12.2813(47) 3135(508) 



TABLE III: Results for the derivative of the phase shift, the derivative of the kinematic function 
(j) P and the Lellouch-Liischer factor for quark masses m/m s = 0.02/0.04. The Lellouch-Liischer 
factor is presented here without the overall V 2 factor (see eq. dU). 

evaluate (fi p '(q*) in the Lellouch-Liischer factor ([[]) at the same point, (c/j + qi_i)/2. 

The results for the derivative of the phase-shift, together with the derivative of the kine- 
matical function <p p which is also an ingredient of the LL factor (see ([I])) and the LL factor 
are presented in Tables Hill and [TV] and plotted in Figures |6] and [7J Our final results for the 
Lellouch-Liischer factor have errors in the range 5-10% when the light quark mass is 0.04 
and 10-20% when it is 0.02. The final error in the LL factor comes predominantly from that 
in 85 /dq*, which contains the finite- volume corrections due to the interactions of the two 
pions. We see from Tables UTTl and HVl that although d<f) p /dq* is larger than 85 /dq* they 
have the opposite sign, so that the effect of the latter is indeed significant. 

V. CONCLUSIONS 

In this paper we have proposed a method for calculating the derivative of the s-wave 
/ = 2 7T7T phase shift directly in lattice calculations. This enables us to evaluate the Lellouch- 
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q* d5/dq* d(j) p /dq* LL factor 



0.20572(73) 


-1.40(26) 


11.3786(67) 


7204(159) 


0.20701(70) 


-1.41(25) 


11.4698(63) 


7095(143) 


0.20825(68) 


-1.43(24) 


11.5585(58) 


6988(137) 


0.20943(66) 


-1.46(22) 


11.6444(55) 


6882(129) 


0.21028(64) 


-1.49(23) 


11.7284(51) 


6817(126) 


0.21083(63) 


-1.50(22) 


11.7684(50) 


6769(116) 


0.21161(62) 


-1.54(22) 


11.8062(48) 


6679(125) 


0.21346(58) 


-1.64(22) 


11.8767(43) 


6448(125) 


0.21580(54) 


-1.87(26) 


12.0750(37) 


6181(152) 


0.21717(51) 


-2.15(34) 


12.2271(35) 


5965(200) 


0.21781(50) 


-2.49(43) 


12.2810(35) 


5737(252) 


0.21833(49) 


-3.86(82) 


12.3196(35) 


4915(476) 



TABLE IV: Results for the derivative of the phase shift, the derivative of the kinematic function 
(j) P and the Lellouch-Liischer factor for quark masses m/m s = 0.04/0.04. The Lellouch-Liischer 
factor is presented here without the overall V 2 factor (see eq. dU). 



Liischer factor which relates the finite- volume K — > tttt matrix elements to the corresponding 
physical decay amplitudes at the quark masses and lattice volumes being simulated. The 
method relies on the use of partially twisted boundary conditions, allowing for the two-pion 
energies, and hence the corresponding phase-shifts, to be evaluated as a function of the 
centre-of-mass relative momenta. 

The feasibility of the method was successfully tested in an exploratory computation. 
The phase shift was calculated as an almost continuous function of q* and the correlations 
between the points in Fig. [5] allows a reliable determination of the LL factor as illustrated 
in figs. [6] and [71 Although the volume of our lattice was small (L ~ 1.9 fm), this study 
does give us confidence that the technique can now be applied to a more physically realistic 
computation of AI = 3/2 K — > tttt matrix elements, such as that currently being undertaken 
by the RBC-UKQCD collaboration (20 ] . 

We see from Tables ILTT1 and HVl that although 9 varies from to it and the total momentum 
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FIG. 6: The derivative of the phase-shift, 4> p ' and the Lellouch-Liischer factor for quark masses 
m = 0.02 and m s = 0.04. a = SkuikE 2 /q* 2 (see eq.([T]) and, as elsewhere in this paper, we do not 
include the overall factor of V 2 in our numerical results). 



varies from 27r/L to 0, the boost back to the centre-of-mass frame results in a very limited 
range of values of q* and hence we only evaluate the phase-shift and its derivative in this 
limited range. Nevertheless, this method allows us to evaluate the Lellouch-Liischer factor 
for AI = 3/2 K — > tttt decays for any kinematics which can be used with periodic or 
antiperiodic boundary conditions. For example, one could imagine calculating the K — > titt 
amplitude with total momentum 2ir/L and then use the technique proposed in this paper 
to determine the corresponding Lellouch-Liischer factor. In this case, the derivative would 
be evaluated from the behaviour of the phase-shift with small twisting angles. As a second 
example imagine introducing anti-periodic boundary conditions for the u quark in the in- 
direction say, and calculating the matrix element between a kaon at rest and two pions with 
momenta ±7r/L. In this case one would determine the derivative from the behaviour of the 



phase- shift with twisting angle 9 Z c 



OSe tO 7T. 



Disappointingly, as explained in 11], the breaking of isospin symmetry induced by choos- 



ing different boundary conditions for the u and d quarks means that the method cannot be 
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FIG. 7: The derivative of the phase-shift, <p p ' and the Lellouch-Liischer factor for quark masses 
m = m s = 0.04. a = 8irV 2 mKE 2 /q* 2 (see eq. (TT]) and, as elsewhere in this paper, we do not include 
the overall factor of V 2 in our numerical results). 

applied to AI = 1/2 transitions. In that case, even in the free theory, the 7r + 7r~ state has a 
different energy from the 7r°7r° state. 
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